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ABSTRACT 

The discovery potential of the Tevatron CDF for the rare B-decay Bg /^^^^ is analysed. 
We find that with an integrated luminosity of 2 fb~^, and using CDF as the example detector, a 
5a combined discovery reach of the Tevatron is possible if the Branching ratio for Bs is 
(1.7 lb 0.46) X 10^''. Such a possible signal for the decay Bg /^^/U" will invite large tan/? values 
and set an upper bound on the heaviest mass of the MSSM Higgs sector in a complete analogy to 
the upper bound of the lightest observable supersymmetric particle set from the excess over the 
SM prediction of the muon anomalous magnetic moment. If for example, the decay Bg 
is found at Tevatron with Br(i?s = 2 x 10^'' then the heaviest Higgs boson mass in the 

MSSM should be less than 790 GeV for tan /3 ^ 50 provided that the CKM matrix is the only 
source for (s) quark flavour changing processes. 
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1 Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) the Higgs bosons are not scalar 
superpartners of known Standard Model (SM) fermions. Apart from a naturally light 
Higgs boson, the rest of the Higgs scalars may well have masses at very high energies 
despite fine tuning matters. We argue here that if the Tevatron observes Bg — * fJ''^IJ'~ , the 
heaviest Higgs boson mass of the MSSM is bounded from above to be less than 1 TeV for 
tan P ^ 50, provided that the Cabibbo-Kobayashi-Maskawa (CKM) matrix [1] is the only 
source for (s) quark flavour changing neutral current (FCNC) processes and CP- violation. 

The pubhshed Tevatron/CDF physics results with luminosity 171 pb""^ is the bound on 
B, ^ [2] 

Bt{Bs //V") < 5.8 X 10"^ , at 90% C.L . (1.1) 

The predicted ratio for Br(S^ in the SM is (3.8 ± 1.0) x 10"^ [3]. The observable 

Bg yU+yU" makes itself special in the MSSM where it is enhanced by orders of magni- 
tude. This enhancement is due completely to the MSSM neutral part of the Higgs sector. 
In particular, Bt{Bs — > IJ'^IJ'^) is proportional to the sixth power of tan/3 and inverse pro- 
portional to the fourth power of the CP-odd (or CP-even) Higgs boson mass. Ma [4]. The 
idea here is that, since the parameter tan/9 is theoretically^ bounded from above and the 
amplitude Bg /i+yU" exhibits a non-decoupling behaviour in the large SUSY mass limit, 
any observation of Bg — ^ fi^fi~ at the Tevatron will set an upper bound on the neutral 
heavy Higgs masses. We present a quantitative analysis of this upper bound in the MSSM 
under the assumptions of conserved R-parity symmetry and no extra source for FCNC's 
other than the CKM matrix itself. 

Although the rare decay Bg —>■ fJ^'^lJ'~ , and the SUSY contribution to the muon anoma- 
lous magnetic moment 5a^, go hand in hand in the minimal Supergravity scenario [17] 
this is not true in a general unconstrained MSSM which we explore here. It is however 
possible, and has been shown in Ref. [30], that following an argument along the lines of 
the previous paragraph but for Sa^, one could set an upper mass bound on the lightest 
observable supersymmetric particle. For complementarity purposes, we repeat the analysis 
of Ref. [30] with the up-to date (and persistently non-zero) measurement of (5a^. 

2 Discovery potential for Bg at Tevatron 

In order for Bg^d ^~ to actually be observable at the Tevatron it would need to ap- 

pear as a clear, unambiguous signal above background events. Because this is an important 
signature, potentially indicating unequivocally the existence of beyond the SM physics, a re- 
quirement that any observed signal must be at least five standard deviations above the level 
of background fluctuations would likely be imposed by the experimenters. Consequently, 
a study of the discovery potential of any future experiment requires that we estimate the 

^ Large tan/3 > 50 — 60 values lead to non-perturbative evolution of the top or bottom Yukawa couplings 
at higher energies. 
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amount of background in the experiment as the data increases. Because a detailed descrip- 
tion of the background levels, and their uncertainties, exists for the CDF collaboration, 
what follows will draw heavily from their observed background levels [2,31,32]. 

To estimate the number of signal events and precision obtained in a hypothetical future 
experiment, we assume that, with no changes in selection criteria for the events from [2], 
the number of background events in the region of interest around the Bs^d masses will be 
of the following form: 

Nbkg = (^bkg ■ X , (2.1) 

where x = J C dt and abkg is the equivalent cross section of the background in the mass 
region of interest. We set the discovery level for the number of Bg^d /^^/^~ candidates 
(Nsig) to be Nsig > ^\/Ni)kg- Using these relations we present the data in the form 

aBs4 ■ Br B,,d ^ = 7^ — ^ y= — , 2.2 

2e ■ a ■ X 2e ■ ai/x 

where e ■ a is the efficiency of the entire selection process multiplied by the detector ac- 
ceptance. The uncertainty in e ■ a is independent of the integrated luminosity so this 
uncertainty has been included as part of the constant error on abkg- ^{x) represents the 
fact that all experiments have uncertainties and that often these uncertainties are based in 
part on the data themselves. As a result we expect the precision on part of the background 
estimate will improve as the data collects. The factor of two accounts for the fact that the 
experiment does not distinguish between the charge conjugate decays. 

The advantage of using equation ()2.2|1 is that all the quantities on the right are obtained 
from a single experimental analysis. In contrast, on the left side of the equation, Br(Bs,d 
fi'^fi~) is the quantity of interest and aBs,d = ^(j{B^) depends on B meson fragmentation 
measurements from other experiments and a cross section measurement in a different decay 
mode (and for this paper, at a different y/s). So equation ()2.2|1 allows the reader to 
more easily adjust estimates of experimental reach for different integrated luminosities, 
fragmentation estimates, and B meson or b quark cross section measurements. 

Henceforth we will focus on the decay Bt{Bs f^^f^~) only. In 171 pb~^ the CDF 
experiment reports Nbkg = 1.05±0.31 events. The largest uncertainty is actually luminosity 
dependent because it comes from the statistical fluctuations of 14 events that remain in the 
unblinded region after some of the final analysis cuts were relaxed. Nbkg can be converted 
into an equivalent background cross section using (Xbkg ± ^(a^) = Nbkg/xo = 6.14 x 10^^ ± 
2.9 X 10"^ ± ^ pb. For this paper xq = 171 pb'^ while e-a = 0.020 ± 0.002. 

Figure ^ shows how the a ■ Br(i?s fJ'^f^~) would scale with increasing integrated 
luminosity given that this decay is sufficiently prolific to allow a discovery using parameters 
from the CDF experiment. The uncertainty band shown here refiects the uncertainty in 
the background estimate and takes into account the additional statistical uncertainty that 
would be added in quadrature if this decay was seen at the 5a level. 

Taking the 5a values from Figure Q for 0.5, 2.0, and 10 fb~^ and using^ j- = 

fx is the probability that a b-quark fragments into a given B meson. The most recent results from the 
Particle Data Group are used, http : / / pdg . Ibl . gov/) 
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Table 1: Tabulated are discovery branching ratios for CDF at three different integrated 
luminosities. The uncertainties are from the background estimate in [2]. 



JCdt, fb-1 


ass ■ Br(i?s /U^/W"), pb 


Br(S, X 10"^ 


0.5 


0.44 ±0.14 


4.8 ± 1.7 


2.0 


0.22 ±0.05 


2.4 ±0.65 


10.0 


0.098 ±0.014 


1.1 ±0.24 



and a{B+) = 3.6 ± 0.6 /ib [33] we obtain the results in Tabled for the branching ratio 
By{Bs —>■ ii^ijT). The analysis above makes use of the detailed reports on this decay 
from the CDF experiment alone. In order to obtain a Tevatron estimate one would need 
to include the D0 data as well. There is a preliminary estimate for an upper limit of 
Br(5s fi^fi') = 4.6 X lO"'^ at 95% CL on D0's public analysis web site [34]. 

It is difficult to account for all of the differences between the two experiments, some 
of which improve and some of which detract from the relative sensitivity between the two 
experiments, but it is clear that a limit on this decay mode similar to CDF's has been 
obtained. Consequently we felt justified in lowering our final estimate of Tevatron reach by 
a factor of -\/2 but present only the CDF-derived values in the paper because of a better 
understanding of that detector by one of the authors. 



3 Bounding the MSSM Higgs and SUSY sectors from 
above 

As we stated in the introduction, excesses of the experimental data over the SM expectation 
for observables like the muon anomalous magnetic moment a^, and the rare B-meson 
decay Bg /i"*"/^", would result in bounding from above the SUSY and the Higgs sectors, 
respectively. In this section, we quantify these statements. 

At large tan/3, the leading supersymmetric contribution to Bt{Bs Z^"*"/^") is depicted 
in Fig.(j2I). Its amplitude depends mainly on 6 beyond the SM parameters, namely, the CP- 
odd Higgs mass (or the CP-even Higgs boson mass because Ma — Mh), the ratio of 
the vacuum expectation values, tan /?, the soft supersymmetry breaking trilinear coupling 
At and the two scalar superpartners of the top quark masses, m^^ and m^^. It is well 
known [4], that the self energy diagram of Fig.Q, enhances the ratio Bt{Bs /W"*"/!") in 
the MSSM by many orders of magnitude. It is worth making a parenthesis here in order 
to understand why this is so. 

The leading contribution in Br(i?s fi^fi^), comes from the amplitude in Fig.(j21). Its 
calculation can be sketched with the following steps 

1 1 

A{bR SLfiRfii) ~ Vb — Vb Vtb fi yt At Vu yt Vt, f{fi, mt^,mtj — j (3.1) 
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Vtfcl^t.l/Ml/^'l/t (^^j (3-2) 



^ V^ffcl//, mf, m^^ , (3.3) 

where \^ is the CKM matrix. The function /(/x, m^^, mj^) is a loop factor which, in the 
hmit of heavy squark masses, behaves hke 1/m?. m? denotes the geometric mean squark 
mass, i.e m? = mi^irti^. This hmit is taken when passing from Eg. 1)3.11) to Eq. ()3.2|) . The 
approximation symbol "~" stands for constant, 0(1), numerical factors irrelevant to our 
point here. 

Eq. fl3.2j) . shows that the amplitude is proportional to one generic (i.e., not originated 
from the conversion of a Yukawa coupling to a quark or lepton mass) tan /3 parameter. The 
other two powers of tan (5 arise in Eq. ()3.3|) from the fact that the Higgs field mediating the 
amplitude is the heavy one and at large tan/3, is [yi,^^ oc mb_^tan/3]. The reader should also 
notice that the light Higgs boson coupling to the bottom quarks or to the leptons is not 
proportional to tan/5 and thus leads to subdominant contributions. It is very interesting 
to note that, in the SUSY decoupling limit, ~ ~ — > cxd in Eq. fl3.2j) . the "Higgs 
penguin" does not decouple from the MSSM. This non-decoupling property makes the rare 
decay Bg —>■ visible even for multi-TeV superparticle spectrum. 

To summarise, squaring the amplitude in ()3.3|) we obtain 

MBs - /iV") oc ^ . (3.4) 

Thus given an upper "theoretical" bound for tan (3 or measuring tan (3 by other means, an 
observation of the decay Bg — ^ /i"*"/^" at Tevatron such as the one we described in Table 1, 
will set an upper bound on Ma- 

In our numerical analysis, we make use of two independent calculations for Br(i?s 
fi~^fi~). The first contains the 1-loop calculation, including corrections after the electroweak 
symmetry breaking [9] with the method of tan/3 resummation as described^ in [17]. Our 
second calculation is based on the effective Lagrangian technique [6,11,13,14]. In particular 
we follow the work [13] where a general resummed effective Lagrangian for Higgs-mediated 
FCNC interactions in the MSSM has been derived, without restrictions to particular as- 
sumptions that rely on the quark- Yukawa structure of the theory or CP-conservation. For 
large tan/5 the two calculations are in good agreement. 

Our calculation of supersymmetric contributions to muon anomalous magnetic moment 
5a^, follows the work in Ref [36]. Recently, two-loop tan /3-enhanced corrections have been 
carried out [37]. These corrections account for ^ 1 x 10~^° if the sparticles are heavier 
than approximately 400 GeV and they are neglected in our analysis here. Finally, for the 
Higgs mass bound we use the numerical code from [38]. 

^Basically, one has to replace the two bottom Yukawa couphngs yt, in Fig.© with the effective yb/{l + 
Arrih) where Ami, cx fitanf] depends on the gluino and sbottom masses [35]. Notice that in Ea. (|3.2() only 
one bottom yukawa coupling survives at the end. 
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We use a high statistics scan over the unconstrained MSSM^ parameter space in the 
region where the soft breaking masses are below 2.5 TeV. We assume that the FCNC cou- 
phngs originate only from the CKM matrix. Furthermore, we assume that the Lightest 
Super symmetric Particle (LSP) is stable and neutral^. That means the LSP is either a 
neutralino or sneutrino. For each set of outputs we define the Lightest Observable Su- 
persymmetric Particle (LOSP) as the second lightest SUSY particle or the third if the 
first two LSP's are both neutralino and sneutrino. Finally, we scan only over the positive 
sign of the /i-parameter consistent with the bound from b — > 57. In any case, results for 
Bt{Bs fJ'^f^ ) ~ 10~^, are independent of this assumption [see Eg. ()3.2|l ]. The unification 
assumption Mi = M2/2 has been assumed for simplicity. We have imposed the current 
experimental bounds from the direct SUSY searches at LEP and Tevatron [39]. For the 
light Higgs boson we have set a rather conservative constraint > 90 GeV since we work 
in the large tan (3 regime. 

Our results are summarised in Fig.Q. Fig.(j2K) displays the ratio Bt{Bs /^^yU^) 
with respect to the CP-odd Higgs mass, M^. Clearly, Br(i?s can be three or 

even four orders of magnitude bigger than the SM expectation and certainly within the 5cr 
Tevatron/CDF reach as shown in the figure, for 2 fb~^ (see Table 1 for other choices of 
luminosity). Even for very large Ma ~ 1 TeV, the ratio can still be enhanced by an order 
of magnitude relative to its SM prediction. However, as Ma increases, Br(i?s — ^ f^~^f^~) can 
not be bigger than certain values. The envelope contour in Fig.(jHK) is well approximated 
with 



tan/?\6 .650 GeV 
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Bt{B, /iV) = 5 X 10-^ (^^j ( ) + 1.0 X 10-« , (3.5) 

We have checked that Eg. ()3.5j) fits well for larger values of tan/3. Now, from Eg. ()2.2j) . if 
CDF with 2 fb~^ (10 fb~^) detects 17(88) signal events for the B, — s> fi^fi that means that 
By{Bs Ai+/i^) 2 X 10^'' and the heaviest Higgs Boson mass will be less than 790 GeV 
for all tan /3 ^ 50. 

A more "natural" prediction for Bt{Bs —>■ fJ'^fi^) is approximated with the inner curve 

= 5 . 10- (^)' (^)^ + 5.0 X 10- , (3.6) 

and represents the statistically populated area of points. Using Fig.(jnti') and our analysis 
shown in Table 1 for 5cr Bs discovery at Tevatron/CDF, we conclude that for 

tan/? < 50, 

Ma < 660 (305) GeV with x = 0.5 fb"^ 
Ma < 790 (360) GeV with x = 2 fb"^ 

Ma < 970 (440) GeV with x = 10 fb"^ , (3.7) 

"^Unconstrained MSSM means here that we impose no restrictions on the boundary conditions for squark 
and slepton masses or trihnear soft SUSY breaking couphngs at the GUT scale. We do not also require 
radiative electroweak symmetry breaking. 

^This is only necessary for the R-parity conserved (we assume here) models. 
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where the numbers in parenthesis fall into the Eq. ()3.6|) . Further uncertainties due to the 
Bg decay constant Jb^, are also reduced when calculating the M^, as compared to the ones 
involved in the calculation of Bt(Bs — > fi^fj.~), since it is sj fss vs. /^^. An error of 5-10% 
on Ma can still be present but can be further reduced once AM^ is measured, possibly at 
Tevatron, along the lines of the last paper in Ref. [3]. Eq. ()3.7|l is the main result of our 
paper. 

The current BNL [40] experimental result on the muon anomalous magnetic moment, 
shows a 2.7(7 difference relative to the SM expectation. The deviation from the SM value 
= a^^P - a™ at 90% C.L reads [41] 

12.9 X 10^^° < < 36.0 x 10"^° , (3.8) 

with a central value Sa^ = 24.5 x 10^^'^ [41] ^. A possible excess can originate from six 
supersymmetric mass parameters, namely, the gaugino masses Mi and M2, the Higgsino 
mass parameter /i, two supersymmetric scalar partner masses of the muons, m^^ and m^^, 
and tan j3. Then, in the mass insertion approximation, one may write this contribution 
simply as [30] 



o2 

Sa^ ~ ml fiM, tan (3 F , (3.9) 



with F oc l/Mg^gY being a loop factor depending on the supersymmertic masses 
(charginos, neutralinos, smuons). Feng and Matchev in [30] used the excess on the muon 
anomalous magnetic moment to set an upper mass bound on the Lightest Observable Su- 
persymmetric Particle (Mlqsp)- We repeat their analysis in Fig. (03), where we plot the 
supersymmetric contribution 5a^, versus the mass of the LOSP. The current excess on 6a^ 
already sets an upper bound Mlqsp ^ 700 GeV, on the LOSP mass. This is clear from 
the envelope solid contour in Fig.(jnb), which is described by the equation 



Sa, = 18 X 10-'" ^ ^^^j . (3.10) 

Thus the current central value [40,41] 6a^ = 24.5 x 10~^° implies that the lightest observable 
Supersymmetric particle should weigh less than 450 GeV. The analogy with the ratio 
By{Bs — >■ /i"*"/^") in Fig.(jnt) is obvious. 

Notice that, from the Fig.(jnK) one cannot set a lower bound on the CP-odd Higgs mass 
from the Tevatron/CDF experimental bound alone. The parameter which can be excluded 
by CDF is the ratio of the soft SUSY breaking trilinear parameter with the CP-odd Higgs 
mass. At /Ma- This bound can be understood easily from Eq. ()3.2|) and can be derived from 
Fig.jSfc) to be 

^ < 12 for tan (3 > 50. (3.11) 

Ma 

Combining Eqs p.llj) and ()3.7|) we arrive at the conclusion that any observation of Br(i?s 
tJi^ IJT) will not only bound the Higgs sector from above but will also set an upper bound to 

^New estimates on the light by hght scattering [42] as well as QED contributions [43] have been 
taken into account when quoting this value. 
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the top trilinear soft SUSY breaking coupling. A soft supersymmetry breaking parameter 
is bounded by a Higgs sector parameter! 

Finally, in Fig.(j21i) we plot the two observables 6a^ and Bt{Bs together. 
Since both are saturated by different mass scales [6a^ by the SUSY and Br(_Bs — > /^^/^~) 
by the Higgs mass scale)], we do not expect any correlation between them in the general 
unconstrained MSSM. This is exactly what we obtain in Fig.(jnii). In fact, the statistical 
distribution of the points with large Br(i?s fJ'~^f^~) and small (5a^ is wider than the other 
way around (we have scanned the parameter space uniformly). However, in the mSUGRA 
scenario, the SUSY mass M^qsp and the CP-odd Higgs mass, M^, are related (mainly) 
through the parameter M1/2. The mSUGRA points in Fig.(jnii) lay along the diagonal in the 
5a^ — 'Bi{Bs fJ'^fJ' ) plane. As it was remarked in Ref. [17], a discovery of Br(i?s — > 
at the Tevatron and a deviation on (as is currently the case [44]) on both observables 
with respect to their SM values will definitely favour some kind of mSUGRA scenario. 

We have checked the robustness of our results presented in Fig.Q against i) a, general 
scan in the multi-TeV region of soft supersymmetry breaking masses, ii) inclusion of various 
supersymmetric phases and Hi) the inclusion from the gluino loop corrections. The last 
two contributions have been implemented in the analysis of Ref. [13]. We find numerically 
that these contributions do not spoil the upper bound of Eq. ()3.7j) . However one expects 
different bounds in the case of non-trivial squark mass intergenerational mixings in the 
supersymmetry breaking sector. The effects of the latter will not be missed at B-factories, 
colliders, or other non-accelerator experiments; where strong constraints on them already 
exist. 

We note in passing, that Tevatron [2] and B-factories [45] search for the Bd leptonic 
decay, Bd ■—>■ /i"*"/^"- Current limits obtained in the region [few x 10~^]. The branching 
ratio Br^Bd — > /i^/i~) is smaller than Bt{Bs — > fi'^fJ^~) by a factor |Vf(i/\4sP ~ 0.06. Any 
observation at Tevatron for Bg IJ''^I^~ is likely to be accompanied with an observation of 
Bd fJ'^f^^ at B-factories if the Higgs sector is responsible for making these leptonic rare 
decays visible. 

4 Epilogue 

We investigate the discovery potential of the Tevatron hadron collider, focusing mainly on 
the CDF experiment, for the rare decay Bg /^"^A*"- We find that Tevatron can reach a 
5a (3(j) observation for 5s-cross section times Bt{Bs f^^f^~) with luminosity given in 
Fig.(^. Based on these results and the knowledge of the cr{Bs) from the Tevatron we show 
in Table 1 the 5a discovery bands for Bt{Bs — > /i^/i~) at several luminosities. For example, 
the Tevatron/CDF with integrated luminosity 2 fb~^, can discover Bs with a ratio 

Bt{Bs At+/U") = (2.4 ±0.65) x 10"^ If this is true, we show that in the MSSM, with the 
assumption that flavour changing processes originate exclusively from the CKM matrix, 
the Higgs sector is bounded from above with the heavy Higgs bosons weight no more than 
790 GeV for tan/3 < 50 [see Fig.(jnK) and Eq. ()3.7|l ]. The situation is completely analogous 
with the muon anomalous magnetic moment searches, where the current excess over the 
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SM prediction seems to indicate an upper bound of 700 GeV on the lightest observable 
supersymmetric particle. 
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CDF Discovery Potential 
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■ Cross section required for 5a discovery 
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Figure 1: Tevatron's CDF discovery potential for the product of the production cross section 
times the Branching ratio [a{pp Bg^d + X) * Bt{Bs fi'^fi~)] vs. integrated Luminosity. 
The red (blue) squares are the center points for 5a {3a) observation. The dashed lines 
delineate the uncertainty in the discovery cross section (both statistical and systematic) as 
a function of integrated luminosity. 
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h,H,A 




Figure 2: The leading supersymmetric contribution to Bg jj,^jj,^at large tan/?. "Crosses" 
and "dots" indicate mass insertions and Yukawa vertices, respectively, and their corre- 
sponding interaction parameters are explicitly given. 
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Figure 3: a) The Br(i?s f^~^f^~) vs. the CP-odd Higgs mass in the MSSM. The lower 
shaded band shows the SM prediction for Bt{Bs ji'^jj,"). The 90% C.L excluded area from 
Tevatron/CDF [2] is also displayed. Following our results from Table 1, the 5a reach at 
2fh-^ forBi^Bs fi^fi ) is also depicted (upper shaded band), b) The SUSY contribution 
to the muon anomalous magnetic moment, 6a^ versus the mass of the Lightest Observed 
Supersymmetric Particle (LOSP). Black points correspond to neutralino or sneutrino LSP. 
The shaded region is the allowed value on at 90% C.L [40,41]- c) The Bt{Bs — > fi^fJ'") 
vs. the ratio of the soft SUSY breaking trilinear parameter with the CP-odd Higgs mass, 
Ai/Ma. d) Combined predictions for 5a^ and Bi{Bs — > fi~^fi~) in the MSSM. For all the 
figures here we have chosen tan/? = 50 and rrit = 178 GeV. 
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